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ABSTRACT 

We present the results of a search for pure rotational molecular hydrogen emission 
from two pre-main-sequence stars, AB Aur and CQ Tau. Observations were made 
using MICHELLE, the mid-IR cchcllc spectrometer at the UK Infrared Telescope. We 
found some evidence for emission in the J = 4 — > 2 line in AB Aur, but no J = 3 — ► 1 
line from either star. We derive upper limits on line flux which are significantly smaller 
than previous line flux estimates based on ISO observations. This suggests that the 
emission detected by ISO is extended on scales of at least 6 arcsec, and does not come 
from the disk as previously thought. 
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1 INTRODUCTION 

Circumstellar disks play a vital part in both star and planet 
formation, and if we are to understand these processes then 
it seems we must understand the composition and physics 
of disks first. The discovery that disks can survive well into 
the main sequence lifetime of a star (after planet forma- 
tion ended in our own system) has ope ned up an excitin g 
window on the planet formation process JZuckermanll200ll) . 
This means that we can directly probe the environment of 
planet formation, even though planets themselves may be 
undetectable. 

The dominant constituent of the disks around pre-main- 
sequence stars is H2. However, observing emission from H2 
which represents the bulk of the warm ( ^200K) gas in th e 
disks is not str aightfoward. Near-in frared (IBarv et all 2002'). 
UV emission (IValenti et alJ 1200 or absorption observa- 
tions jRoberee et al.ll200ll) trace only the hottest gas or 
are dependent on the line of sight through the cirumstel- 
lar disk. Mid-infrared observations of the pure rotational 
emission lines of H2 trace the bulk of the warm gas, but 
are difficult from ground-based sites, due to telluric absorp- 
tion and the high mid-infrared background emission. Thus, 
many previous studies of the gas component of such disks 
has relied on observations of tracer molecules, part icularly 
CO dMannines fc SargentJlioOul iDuvert et alJl2000l) . While 
these observations have been successful in tracing the veloc- 
ity structure in disks, it seems that CO is a poor indicator 
of the total gas content. This is because CO can be de- 
stroyed by UV photons or frozen onto dust grains, leading 
to a severe underestimate of the total amount of gas. Most 
estimates of disk mass are based on sub-mm observations of 
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the d ust continuum feeckwith et alJll99(J : ISvlvester et alJ 
l200ll) . but these rely on assumptions about the composi- 
tion and size distribution of the dust grains, and also on 
the gas/dust mass ratio. More importantly, we know that 
the dust and gas do not evolve in the same way - disks 
around main-sequence stars are composed mainly of dust, 
while disks around pre-main-sequence stars are composed 
mainly of gas. It is important to be able to measure gas 
properties independently from the dust properties using di- 
rect observations of H2 . 

Recently, the Infrared Space Observatory (ISO) has 
been used to search for H2 pure-rotational emission 
lines from pre-main-sequence stars, and also yo ung main- 
sequenc e stars which are known to have dust disks. Frhi et alJ 
fcOOlbl) report detections of the H2 J = 3 — > 1 transition 
(A = 17.035 urn) and H 2 J = 2 -» transition (A = 28.221 
pm) . H2 masses derived from these observations imply large 
masses of H2 in the disks around main-sequence targets 
<Thi et al.ll2001alh reversing the previous belief that dust 
is the main component. However, as the ISO beam is large 
(14 x 27 arcsec) compared with the disks and the observa- 
tions have moderate spectral resolution (R ~ 2000), these 
results cannot distinguish circumstellar disk emission from 
that of a foreground cloud or extended source. 

Follow up observations of both the pure rotational emis- 
sion lines and H2 in absorption have provided strong ev- 
idence that the warm H2 is extended in sources selected 
from the ISO studies. Using the the echelle spectrometer, 
TEXES, iRichter et al.l |2002) obtained observations of six 
sources, three of which have reported ISO H2 detections 
<Thi et al.ll2001bl) : HD163296, AB Aur and GG Tau. They 
obtained upper limits on the line fluxes from five sources, 
with a possible 2 sigma detection of the H2 J = 3 — > 1 S(0) 
line in AB Aur and conclude that the majority of the H2 
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emission detection by ISO is extended on scales of 5 arcsec 
or more. The far-UV spectrum of B eta Pictoris presented 
bv lLecavelier des Etangs et alJ (j200lfl shows a complete ab- 
sence of H2 absorption lines, contrary to expectations if large 
quantities of H2 are present within the edge-on dust disk. 

We present the results of our attempts to confirm 
the ISO results using MICHELLE, the mid-infrared echelle 
spectrometer on the ground-based United Kingdom Infrared 
Telescope (UKIRT). We targeted two pre-main-sequence 
stars, CQ Tau (distance 100 pc) and AB Aur (distance 142 
pc), both of which have good ISO detections and one o f 
which (AB Aur) was also observed bv lRichter et al.l $2002 ). 
In both cases, the dust disk is expected to be unresolved 
spatially by MICHELLE. 



2 OBSERVATIONS AND DATA REDUCTION 

Mid-IR echelle spectroscopy of AB Aur and CQ Tau was 
carried out on 2001 December 29 and 30 using MICHELLE 
JClasse et al.lll997l) on UKIRT. AB Aur and CQ Tau were 
observed in the H2 J = 3 — > 1 transition (A = 17.035 u.m), 
providing direct comparison with the ISO observations. AB 
Aur was also observed in the H2 J — 4 — > 2 transition 
(A = 12.279 \im). The H 2 J = 2 -> transition (A = 28.221 
fun) was used in the ISO observations, but is not readily 
accessible from the ground due to the structure of the at- 
mospheric transmission at 28 urn. Details of the observations 
are shown in Table 1. 

MICHELLE is a mid-IR imager and spectrometer which 
uses a 320 x 240-pixel Si: As array and operates between 8 
and 25 microns. For these observations the echelle grating 
was used, to provide maximum sensitivity to the H2 lines 
which are expected to be unresolved. AB Aur was observed 
with a 2-pixel slit width; CQ Tau with a 4-pixel slit width. 
The resolving power was R — 15 200 at 12 u.m and R = 
18 600 (R = 9 300) at 17 urn for a two (four) pixel wide 
slit; as a result the H2 lines are well separated from the 
bright sky lines. A single order from the grating is selected 
using blocking filters. The pixel field of view varies with 
wavelength due to anamorphic demagnification. The values 
are given in Table 1. 

The point-like sources were all observed by nodding 
along the 70 arcsec long slit, which was orientated north- 
south on the sky. The separation of the two beams was 15 
arcsec and the dwell-time in each beam was 40 seconds. 

Data reduction was done using software written in the 
idl programing language. Pairs of frames from different nod 
positions were subtracted, giving a set of difference frames 
each containing a positive and a negative version of the 
source. In each frame, the positive and negative versions 
were combined to give a single 2D spectrum. Before the set 
of 2D spectra were combined into a final image, the dataset 
was despiked using an iterative sigma cut. To do this, all of 
the 2D spectra were compared, and datum were flagged as 
bad if they were further than a specified number of standard 
deviations from the mean of their pixel. This sigma cut was 
done three times with threshold settings of 6, 4 and 3 sigma. 

Once a final 2D spectrum had been produced, ID spec- 
tra were extracted in two ways: a wide extraction which 
summed all the emission in a 6.5 arcsec region centered on 
the source, and an optimal extraction of the continuum. The 
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Figure 1. Spectra of the H2 J = 3 — > 1 transition for AB Aur 
and CQ Tau. For comparison, gaussians are shown wit h the same 
integrated flux as the ISO measurements reported in iThi et alJ 
assuming an intrinsic line width of 20 km s . The 
CQ Tau spectrum has been masked out at regions of high noise 
(shown by the use of dotted lines); the lower panel shows the 
standard star observation, scaled from 0.0 to 1.0, to demonstrate 
that the masked regions correspond to atmospheric absorption 
lines. For AB Aur, the standard star observation was saturated, 
which causes the variability in the continuum level. 

optimal extraction method takes a weighted mean of the 
pixels around the source, giving more weight where the con- 
tinuum is brighter. This method gives the lowest possible 
noise on the continuum level, but if the H2 emission region 
is more extended than the continuum emission then it will 
underestimate the H2 line flux. For both extraction methods 
a sky subtraction was done first, estimating the residual sky 
signal in each column using pixels between 3.5 and 7 arcsec 
from the source. 

Standard stars were observed and reduced in the the 
same way as the science targets. The final spectrum for 
each object was divided through by a standard spectrum 
to remove wavelength dependant variations in atmospheric 
and instrumental transmission. The standard stars were ob- 
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Object 


Date 


Wavelength 


Slit width 


Pixel scale on-sky 


Integration time 








pixels 


arcsec 


seconds 


AB Aur 


29 Dec 2001 


12.279 


2 


0.59x0.24 


1900 




29 Dec 2001 


17.035 


2 


0.50x0.28 


3400 


CQ Tau 


30 Dec 2001 


17.035 


4 


0.50x0.28 


5900 



Table 1. Summary of observations. The pixel scale on the sky is subject to anamorphic magnification due to the quasi-Littrow mode of 
operation of the MICHELLE spectrometer. The long axis is along the slit. 



Object 


Wavelength 


Line flux 






\Jxn 


xl0 — 14 ergs s _1 cm -2 








Wide extraction Optimal extraction 


ISO Result 


AB Aur 


12.297 


5.4 ± 3.8 4.8 ± 2.1 






17.035 


< 14.6 < 8.9 


30 ± 9 


CQ Tau 


17.035 


< 13.0 < 6.7 


40 ± 12 



Table 2. Summary of results. All upper limits are 3 sigma. ISO results are taken from lThi et alj |2001b). 

served using the same rows on the detector as the sources, 
and division by the standard effectively flat-fields the detec- 
tor. For the 17 [im observations of AB Aur the standard star 
observation was saturated, so we were unable to make this 
correction. This causes variation in the continuum level, but 
would not prevent us from detecting a narrow emission line. 

Wavelength calibration was obtained by measuring the 
positions of sky lines, and comparing the positions to a 
model sky spectrum. We estimate a 1 sigma error on the 
wavelength calibration is about one pixel (~ 3 x 1CP 4 ]xm), 
and is largely due to the difficulty in centroiding the broad 
sky lines. The spectra were shifted to the object's rest frame, 
accounting for the radial velocity of the earth and object's 
heliocentric radial velocit y (+8.9 km s~ for AB Aur an d 
+22.6 km s" 1 for CQ Tau (iBarbier-Brossat fc Figonll999|) 'l. 
Flux calibration was based on the observed conti nuum level 
and pu blished continuum flux measure ments from|TJii_et_alJ 
(2001b) for both the 17 urn spectra and lRichter et alJ 
for the AB Aur 12 |i.m spectrum. 



3 RESULTS 

Figures 1 and 2 show the final spectra produced using an 
optimal extraction. At the location of sky lines the noise 
increases dramatically, so these regions have been masked 
out (i.e. they were not used when estimating the continuum 
level and the overall noise) . On the 17 urn figures a gauss ian 
with the total line flux measured by iThi et al ] fcOOlbl) is 
shown, to allow comparison between the ISO results and 
our spectra. As the ISO observations do not resolve the line, 
we have assumed an intrinsic line width of 20 km s _1 (see 
Section 0] for discussion of this value), and an instrumental 
width of 16 km s _1 for AB Aur and 32 km s" 1 for CQ Tau, 
giving a total width of 26 km s _1 and 38 km s _1 respectivley. 
The MICHELLE sensitvities should have produced a clear 
detection of the line, if all of the flux in the 14 x 27 arcsec 
ISO beam were contained within our extracted beams, as 
expected for emission from an unresolved disk. 

To estimate the line flux, all of the flux within a 
40 km s _1 wide bin centered on the systematic velocity was 
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Figure 2. Spectra of the H2 J = 4 — » 2 transition for AB Aur. 
The vertical dotted line shows to the expected location of the 
transition. Noisy regions have been masked out and the standard 
star is shown in the lower panel, as in Figure 1. 



summed, after subtracting the continuum level. Fluxes and 
upper limits for all extraction methods are shown in Table 
2. There is no evidence of H2 emission in either of the 17 [im 
spectra. The 12 (j.m optimally extracted spectrum gives flux 
of 4.8 ±2.1 ergs s -1 cm -2 , but as this is only a 2.3 sigma 
result should not be regarded as a definite detection. How- 
ever, this resu l t is co nsistent with the tentative detection by 
iRichter et ail J2002t . who measure (2.0 ± 1.0) x 10" 14 ergs 
s _1 cm~ 2 with a width of 10.5 km s _1 . 



4 DISCUSSION 

Despite the enhanced sensitivity to narrow emission lines 
from a point-like source, these new spectra of AB Aur and 
CQ Tau show no convincing detections of either the 0-0 S(l) 
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or 0-0 S(2) lines. In order for the 0-0 S(l) H2 emission to 
be detectable by ISO but not by MICHELLE, the emis- 
sion would have to be either spectrally broad or spatially 
extended. We consider these two explanations in turn. 

The spectral resolving power of the ISO ob s ervatio ns is 
R (= A.) = 2400, and the H 2 lines iThi etaD l)2001bh ob- 
served are unresolved. In calculating our upper limits, we 
assume that the lines had an intrinsic width of 20 km s . 
The dominant line broadening mechanisms for gas in a disk 
are thermal broadening, and bulk motion due to the disk's 
rotation. Ass uming the g as has a temperature of 200 K (as 
estimated by IThi et al.l pOOlbTO . the thermal broadening 
would only be 1.5 km s~ . The amount of rotational broad- 
ening depends on the distance from the star, and the in- 
clination of the disk. If we assume a distance of 2AU (a 
reasonable choice, given the gas temperature of 200 K), this 
gives a width of 15 km s _1 , assuming an inclination of ~ 45° . 
Therefore, 20 km s _1 is a conservative estimate of the line 
width. In fact, the line FWHM would have to be over 100 
km s _1 to make the ISO and MICHELLE measurements 
consistent. 

The possibility of extended emission, due to excitation 
of a remnant cloud of molecular gas by UV photons or by 
shocks, is plausible. Assuming that the ISO measured flux 
is uniformly distributed over the 14 x 27 arcsec aperture, 
the fluxes arc consist ent wit h UV-e xcited emission strengths 
predicted by iBurton et, alJ |[l992) for moderate UV field 
strengths (10 3 — 10 4 times that of the interstellar medium) 
and moderate densities (iih 2 ~ 10 3 — 10 4 cm~ 3 ). Over the 
MICHELLE beam used to obtain the spectra shown in Fig- 
ure 1, the line flux would be an order of magnitude fainter 
than the sensitivity limit obtained. Moreover, the sky sub- 
traction region (3.5-7 arcsec radius) means that our obser- 
vations are insensitive to emission at radius >3.5 arcsec. 
The hypothesis of extended emission is strengthened for AB 
Aur since a map of 12 CO 3-2 shows an extended molecu- 
lar cloud around the centr al source, with a size of at least 
30 arcsec <Thi et al.ll2001bh . For CQ Tau, no 12 CO 3-2 line 
was detected, but this could be reconciled with extended H2 
emission if the CO has been destroyed by photodissociation 
jDent et al.lll995l) . 

In interpreting our results, an important question con- 
cerns the optical depth of the dust disk. If the dust is opti- 
cally thick at the wavelengths we observed and the gas and 
dust are at the same temperature, then the gas would be 
equally likely to absorb or emit a photon, so no emission line 
would be produced. Line emission will only occur if there are 
gaps in the dust disk, or if there is a hot atmosphere above 
the disk. Though these situations are quite plausible, only 
a small fraction of the gas in the disk would be observed in 
this case, and any calculation of the gas mass would under- 
estimate the true mass. It is therefore essential that we know 
the optical depth of the disk in order to correctly interpret 
our results. 

The optical depth of the disks can be estimated from 
radia tive transfer models based on the SED of the disks. 
iNatta et all lj200lTl made models of the disk SEDs for both 
AB Aur and CQ Tau, and we can use their fitted param- 
eters to determine whether the disks are optically thick at 
12 and 17 urn. Based on their estimates of disk mass and 
disk radius, and using the usual assumption that the dust 
opacity at 1.3 mm is /ti.3 mm = 0.01 cm 2 g _1 , the average 
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optical depth of the disks at 1.3 mm is 0.03 for AB Aur 
and 0.24 for CQ Tau. Then, using a power law dust opacity 
index (k oc A - ' 3 , where /3 was estimated to b e 2.0 for AB 
Aur and 1.0 for CQ Tau in INatta et"ai]l200ll) . the optical 
depth can be estimated at shorter wavelengths. This indi- 
cates that the disks become optically thick at wavelengths 
shorter than 200 urn for AB Aur and 300 urn for CQ Tau. 
Both disks must therefore have a high optical depth at 12 
and 17 [im. 

This result has two important implications. Firstly, this 
means that we cannot use our upper limits on line flux to 
calculate upper limits on the mass of H2 in the disk, as most 
of the gas in the disk is hidden by the optically thick dust. 
The H2 line at 12 \im (if confirmed) must arrise from a hot 
atmosphere above the disk, or from gaps in the dust disk. 
Secondly, if the the lines detected by ISO do arise in the 
disk, then only a small fraction of the gas was observed and 
the total gas mass must be much higher than was estimated. 

In conclusion, we have searched for pure rotational H2 
emission lines from the circumstellar disks of AB Aur and 
CQ Tau. We tentatively detect emission from AB Aur, but 
find no evidence for emission from CQ Tau. The likely rea- 
son for this non-detection is that the disks are optically thick 
in the mid-IR. Our upper limits on line flux are significantly 
smaller than previous line flux estimates based on ISO ob- 
servations, which suggests that the emission detected by ISO 
is extended on scales of at least 6 arcsec (or > 60AU). 

Mapping the emission over the area covered of the ISO 
beam is one approach to confirming these conclusions. De- 
tecting the faint extended emission will require a consid- 
erable gain in sensitivity over these observations, such as 
is anticipated using MICHELLE on the Gemini North 8-m 
telescope. 

The United Kingdom Infrared Telescope is operated by 
the Joint Astronomy Centre on behalf of the UK Particle 
Physics and Astronomy Research Council. We wish to thank 
the staff at UKIRT for enabling these observations. We also 
thank Mark Wyatt and Alistair Glasse for helpful discus- 
sions of this work. 
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